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The techniques of shortcuts to adiabaticity have been proposed to accelerate the "slow" adiabatic processes 
in various quantum systems with the applications in quantum information processing. In this paper, we 
study the counter-diabatic driving for fast adiabatic spin manipulation in a two-electron double quantum 
dot by designing time-dependent electric fields in the presence of spin-orbit coupling. To simplify 
implementation and find an alternative shortcut, we further transform the Hamiltonian in term of Lie 
algebra, which allows one to use a single Cartesian component of electric fields. In addition, the relation 
between energy and time is quantified to show the lower bound for the operation time when the maximum 
amplitude of electric fields is given. Finally, the fidelity is discussed with respect to noise and systematic 
errors, which demonstrates that the decoherence effect induced by stochastic environment can be avoided in 
speeded-up adiabatic control. 

ELECTRON spins in quantum dots (QDs)'"'' have been extensively investigated for potential applications in 
quantum information processing, as spins in QDs are expected as a possible realization of qubit in quantum 
information science and technology'. Especially, a two-electron double QD can be further regarded as the 
smallest network to implement quantum computation, in which the highly entangled spin state, i.e. the singlet, 
can be generated. Requirements of precisely controlled qubits have intensively stimulated the detailed studies of 
the interactions in double-dot systems'** and the observations of phenomena thereby, such as Pauli spin blockade" 
and Coulomb blockade'. Furthermore, the demands for achieving efficient quantum computations and avoiding 
decoherence motivate us to manipulate spin states in double QDs in a fast and robust way. There are several 
methods to manipulate spin in QDs, such as electron spin resonance induced by magnetic field oscillating at the 
Zeeman transition frequency' and electric control with spin-orbit (SQ) coupling"'. Recently, conventional "rapid" 
adiabatic passages in quantum optics, for example, Landau-Zener scheme, have been extensively used to spin 
control in single QD'", coupled double QD", tripled QD'^, which can be applied to prepare entanglement states" 
and quantum logical gates, such as NOT'"* and CNOT'^ gates. 

Shortcuts to adiabaticity""" have been proposed to speed up the adiabatic process without final excitation with 
many applications in atomic, molecular, optical physics, many-body physics, and even spintronics, see recent 
review'". In a single QD, we applied the inverse engineering method" to design a fast and robust protocol of spin 
flip in the nanosecond timescale'^", based on the Lewis-Riesenfeld invariant theory'^'. Furthermore, in a two- 
electron QD, more freedom in the applied electric fields provides the flexibility to control spin states by the 
invariant dynamics and controllable Lewis-Riesenfeld phases^^. An alternative shortcut is provided by counter- 
diabatic control proposed by Demirplak and Rice^"*, equivalent to tansitionless quantum driving"'"'. This technique 
was originally utilized to fast adiabatic control in two-level quantum systems theoretically"'"'' "^'' and experiment- 
^jjy25,26 5j^Qj-( afterwards, it has been extended to multi-level systems"'^', and even many-body systems^" 

In this Report, we propose a fast and reliable protocol to generate the entangled spin states by using counter- 
diabatic driving. The external electric fields are designed for rapid spin control in a two-electron double QD in the 
presence of a static magnetic field and SO coupling. We apply the electric fields, instead of magnetic fields, and 
take advantage of SO coupling, since the time-dependent electric fields are easy to be generated on the nanoscale 
by adding local electrodes'. In addition, as comparing to a single QD, counter-diabatic driving is applicable in a 
two-electron double QD, as there exists more freedom with four controllable parameters, x and y components of 
the external electric fields for each dot. To simplify the experimental setup and reduce the device-dependent noise, 
we further apply the concept of multiple Schrodinger pictures'^ to find an alternative shortcut with only x 
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Figure 1 | Schematic diagram of a two-electron double quantum dot in 
the presence of external electric fields and spin-orbit coupling, where the 
singlet state and the lowest one of triplet states are considered as effective 
two-level system, when |/ + A|<§;/ with Zeeman term A = g/isB. 

component of the applied electric fields. Moreover, we also quantify 
how the electric fields increase with shortening the time, to provide 
the lower bound of operation time for a given maximal amplitude of 
electric fields. Finally, the stability of designed shortcuts are discussed 
with respect to decoherence and systematic errors. Our approach 
presents a simple way to manipulate the singlet- triplet transition, 
which could be useful for rapid entanglement state preparation. 

Results 

Two electrons are confined in a double QD, described as a quartic 
potential in Fig. 1, where they are isolated by Coulomb blockade'^. In 
the presence of the applied magnetic fields, the lowest four eigen- 
states of the system can be expressed by singlet and triplet for S = 0 
and S = 1 in the basis of jS, S^. This report presents a method to 
achieve fast adiabatic transition between the triplet and the singlet. 
We design the electric fields mx — y plane to manipulate spin states 
with static magnetic fields along z direction in each dot, considering 
structure-related Rashba (a) and bulk-originated Dresselhaus (/?) for 
[110] growth axis. If the energy difference between the singlet and the 
lowest one of the triplet is much less than the gap between the singlet 
and the triplet, we focus on the state transition between these lowest 
two, as shown in Fig. 1, where Lande factor ^ < 0 like in GaAs and 
InAs QDs. 

By choosing j 1) = (1, 0)'^ and | — 1) = (0, l)'^, referring to the states 
|0, 0) and |l, 1), respectively, we may first take the reference 
Hamiltonian as 



h fZ 



iY 



-iY 



(1) 



where Y=- V2ae{Al-Al) j hc,Z={-J-A)/h+eP{Al +Al) /he, 
and AJ is determined by the electric fields, £j{t) = —{l/c)dAj/ 8t. 



The subscriptions j = L, R represent the left and the right dots, 
respectively. Here we assume the ansatz of the vector potentials is 
AJ=Ao{tanh[[t-ajtf)/[wjtf)] -|- l}, where ai = 0.54,0;; = 0.48, 
Wi = = O.l. The ansatz of vector potentials satisfies the condition 
AJ{0) ~ 0 and guarantees that the electric fields £j start to be driven 
from f = 0, that is, S'' = 0, when f s 0. When the adiabatic condition 
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is fulfilled, the spin state will evolves from | — 1) to 1 1) adiabatically 
along one of instantaneous eigenstates. When the final time is tf = 
1 1 ns, the spin state is completely inverted, and the final population 
of 1 1) is larger than 0.9999. 

Shortening the manipulation time to tf= 2 ns, shrinking Aj into 
this time duration and keeping the same amplitude, we can find the 
state evolution is no longer adiabatic and the final state cannot reach 
1 1) at the final time. The same profiles of time-dependent Y and Z 
terms in Hq are shown in Fig. 2 (a) for different operation times, tf. 

Counter-diabatic driving, equivalent to transitionless quantum 
driving"''^'' "''', provides supplementary time-dependent interactions 
Hi to cancel the diabatic couplings of Hg, and make the process fast 
and adiabatic, where Hi is" 



h fO X 



2\X 0 



(3) 



with X = \/2ae (A-[ —A^g^j he, driven by if the difference between y 

component of two electric fields. As a result, the exact dynamical 
evolution of total Hamiltonian H = Hg + coincides with adiabatic 
approximation of the reference Hamiltonian Ho- However, to imple- 
ment accelerated adiabatic transitions more energy price has to pay, 
that is, the maximal amplitude of in the X term increases when the 
finally time tf is shortened. This can be intuitively understood from 
time-energy uncertainty principle, that is, AJ is proportional to 1/tf. 
Since £j{t) = —{l/c) 8 A, /dt, the larger value of £J and if are finally 
required for the shorter time, tf, as shown in Fig. 2 (c). 

In reality, the electron spin is subject to the device-dependent 
noise, which could be the amplitude noise of the electric fields^". It 
can be quite important, especially when the electric fields are rela- 
tively weak. From the above analysis, we find that four controllable 
parameters, £J and x and y components of the electric fields 
for each electron in a double QD should be applied. If y compo- 
nent of the electric fields can be reduced, we can remove the ampli- 
tude noise from y component of the electric field. In addition to 
decreasing the total decoherent effects resulting from the device- 
dependent noise, the cancellation oiy component of the electric field 
might be also useful to simplify the setup. To this end, we apply 
the concept of multiple Schrodinger pictures to find an alternative 
way to implement the shortcuts. Making unitary transformation of 
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Figure 2 | (a) Time dependence of 7 (solid blue line) and 2 (dashed red line) terms of Hq. (b) The applied electric fields E\ (solid blue line) and E\ (dashed 
red line) drive the state transition of Hq adiabatically, with tf = 11 ns. (c) The applied electric fields (solid blue line), i;^ (dashed red line) andfj^, (dot- 
dashed green line) drive the state transition of H in a fast adiabatic way with shorter time tf= 2 ns. 
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Hamiltonian H^^-" by a rotation around z axis with the angle n/2 ■ 
we obtain 



H = 



+ (j> iQ 
-iQ -Z- 



(4) 



without term, where tan (f) = Y/X and Q = \/X^ + Y^. Again, the 
maximal amplitude of Q will increase when decreasing time tf, due to 
the fact that X becomes dominant (the maximal amplitude of Y is 
unchanged). The Hamiltonian H is equal to the original one Halt = 
0 and tf, which guarantees that the initial (final) states of H and H 
coincide. However, the dynamics is not same during the intermediate 
process, although the populations are always equal. Accordingly, we 
may acquire two new controllable x component of the electric fields, 
£1" and fj", calculated from Eq. (4), see Fig. 3. 

Discussion 

Comparisons of £^ and £^ provided by different times suggest that 
stronger electric fields have to be used for shorter times, though the 
amplitude of electric fields might be optimized by using superadia- 
batic iterations''"'. However, the amplitude of electric fields cannot be 
arbitrarily large simply because strong fields may destroy the sys- 
tems. In order to quantify the energy price mentioned above, we 
demonstrate the relation between the maximal values of electric 
fields and the operation time tf, see Fig. 4. The maximal amplitude 



of electric fields, £^.^^ = max[\£'^ 
very short times, 

£max ^ 



\£^"\), fulfills the scaling law at 



(5) 



since £j" oc AJ ^ tf and ccl/tf go to infinity in the limit of tf^ 0. 
The asymptotic exponent of implies that the minimal time should 
be oc £^ll^, which provides the lower bound of operation time when 
the maximal amplitude of electric fields is given. If the spin system in 
quantum dot, rather than the atom in harmonic trap, is considered as 
working medium in the cooling cycles of quantum refrigerator, the 
minimal time for the (accelerated) adiabatic process, bounded by 
the energy, could be relevant to the third law of thermodynamics 
and the unattainabUity principle" '^. 

For a realistic setup, the coupling to the stochastic environment is 
a general scenario to be considered, where the hyperfine interactions 
with the nuclear spin could play important role at low temperature. 
To study the decoherence effect, we present the master equation for 
the density matrix"" in a generic form: 



P= ■ 



HA 
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Figure 3 | Electric fields of (soUd blue line) and E'j^ (dashed red line), 
designed from the Hamiltonian H, see Eq. (4). 
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Figure 4 | Dependence of E^-^^y^ on short time tf (solid blue line), where the 
dashed straight line shows the asymptotic exponent of ty; i.e. fmaxOC 1 jif. 

where y is the dephasing rate. Solving the Bloch equation, we can 
obtain the final fidelity (P = pn) for different times, see Fig. 5, and 
demonstrate that the faster manipulation increases the fidelity with 
less influences attributed by decoherence. 

To demonstrate the feasibility of our protocol, we also check the 
stability with respect to systematic errors in ipj". The real electric 

fields can be £-f''' = £-f(l + l), where / is the relative deviation. 
The dependence of fidelity _F on a is exhibited in Fig. 6 for different 
times. Different from decoherence affected by the stochastic envir- 
onment, fidelity is more stable with larger tp since the systematic 
error considered here depends on the amplitude of electric fields. 
In general, the speeded-up adiabatic protocol has different stability 
with respect to different types of noise and systematic errors. 
Alternatively, one can combine the inverse engineering and optimal 
control theory to pick up the most robust protocol in quantum two- 
level systems in presence of different noise and errors"". 

Methods 

Effective Hamiltonian. The total spin-dependent Hamiltonian consists of 
Heisenberg term, Zeeman term, and interactions between the electric fields and the 
electrons, expressed as 



(7) 



The subscripts j — L,R represent the left dot and the right one, respectively. Zeeman 
term is A — gfiBB with the equal magnetic fields B applied to the left dot and the right 
one in z direction, and A, are the vector potentials of the electric fields. The spin 
operators of two electrons confined in each dot are s, — (T,/2 with z component 5J. The 

Heisenberg term Jsi • Sn describes the exchange coupling / between two spins. The 
example of a double QD of GaAs-based structure {g—— 0.44) is taken with B = 3.7 T. 
The energy gap between the singlet and the triplet is / — 0.1 meV, so that 
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Figure 5 | Fidelity F versus dephasing rate y with respect to = 2 ns 
(solid blue line), f/ = 3 ns (dashed red line), = 4 ns (dot-dashed black 
line). 
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Figure 6 | Fidelity inversus /i with respect to = 2 ns (solid blue line), t/= 
3 ns (dashed red line), tf = 4 ns (dot-dashed black line). 



|/-t-A|// = 0.06<Cl. SO coupling term of Hamiltonian includes structure-related 
Rashba (a) term and bulk-originated Dresselhaus (/?) term for [110] growth axis, 

j J 

SO that the spin- dependent velocity operators become 

v;'"=|k„..(,)i. 



(8) 



(9) 



As a result, after shifting some quantity of Htotah we can derive a 2 X 2 Hamiltonian 

(10) 

where Z, 7 are A|-dependent while X is -dependent, seen in the section above. 



hf Z X+iY 
2 \X-iY -Z 



Counter-diabatic driving and Z-axis rotation. Naturally, we separate the 
Hamiltonian H into two parts, Hq and Hi , where Hq includes the Y and Z terms driven 
by the x components of electric fields applied in each dot, and Hj includes only X term 
driven by the y components. The strategy of counter-diabatic driving in a two- 
electron double QD is to set Hq as reference first, which could be not adiabatic at all. 
Next, we calculate and add the complementary interaction Hi to cancel the diabatic 
couplings of Ho and make the spin control fast and adiabatic Actually, the 
separation of Hamilton H (10) into Hq and Hj depends strongly on the choice of 
growth axis [110]. For instance, if the growth axis [111] is chosen, the SO coupling 
term should be modified as 



and the 2X2 Hamiltonian (10) becomes 

h f -{J + A)/h X + iY \ 



H- 



X-iY 



{1 + ^)1^)^ 



(11) 



(12) 



with X- \/i'y.e{A{ -A^^) j he and - \/2(a + /?)e(A^ -A^) /he. Therefore, the 

approach presented here is not valid, since the reference Hq and the counter-diabatic 
driving Hi can not be naturally separated and calculated. 

Here counter-diabatic driving is applicable in a two-electron double QD, as in 
Hamiltonian H (10) there exists freedom with four controllable parameters, x and y 
components of the external electric fields for each dot. This is different from the 
Hamiltonian in a single QD where there are only two controllable parameters, x andy 
components of the electric field, so that it is impossible to produce the required aU- 
electrical interaction by counter-diabatic driving^". 

Furthermore, using multiple Schrodinger pictures to describe various physical 
settings sharing the same dynamics is helpful to find alternative shortcuts, when the 
counter-diabatic term is difficult or impossible to implement". One can transform 
the Hamiltonian based on Lie algebra to cancel the unwanted component of 
Hamiltonian^". Applying this concept, we make unitary transformation of 
Hamiltonian H by z-axis rotation. While the original dynamics satisfies 
ihdi^{t)^im'{t), the new dynamics is given by ihdt^'{t)^H^'{t), where 
^'{t)^U^'^{t),H^U'{H-K)U^ndK^ihiJU\ In our case, we use the unitary 
operator U = ihl2)e~ ''^''^^~'^^''% to obtain the Hamiltonian (4). 
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